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In the case of massless current quarks we find that the 
breaking of chiral symmetry usually triggers the generation 
of an anomalous magnetic moment for the quarks. We show 
that the kernel of the Ward identity for the vector vertex 
yields an important contribution. We compute the anoma- 
lous magnetic moment in several quark models. The results 
show that it is hard to escape a measurable anomalous mag- 
netic moment for the quarks in the case of spontaneous chiral 
symmetry breaking. 



Theoretically, the various hadronic electromagnetic 
form factors are usually described in terms of pole dom- 
inance together with contributions arising from virtual 
mesonic exchanges |l]] . A third contribution to the elec- 
tromagnetic form factors should come from the quark mi- 
croscopic interaction itself, in close analogy with QED. It 
is clear that these three scenarios should not be indepen- 
dent but just three different aspects of the same model. 
This desideratum can be achieved, at least qualitatively, 
in terms of a quark field theory displaying spontaneous 
breaking of chiral symmetry, (SxSB). In such a descrip- 
tion any hadron, when seen from the trivial vacuum Fock 
space, appears as a collection of an infinite number of 
quark antiquark pairs together with the appropriate va- 
lence quarks. It happens that the contributions of this 
quark sea can be summarized in terms of a new set of va- 
lence quasiquarks which now carry the information on the 
details of the physical vacuum trough a modified prop- 
agator p|. In this fashion we recover the simplicity of 
the constituent quark picture. It is the role of the Ward 
identities to ensure charge conservation throughout this 
process. And this they do at the expenses of the quark 
magnetic moment which, in general, becomes non-zero. 
As will be shown in this paper, to maintain throughout 
the process of SxSB a zero anomalous magnetic moment 
for the quarks, constitutes the exception rather than the 
rule and is just the consequence of particular choices for 
the Lagrangian. However, the B.C.S. diagonalization of 
the Hamiltonian, (mass gap equation) does not preclude 
quark pair creation or anihilation processes to occur. In 
fact it sets the strength of mesonic contributions for such 
physical processes as decay widths and meson- nucleon in- 
teractions ||, among others. The counterparts of these 
processes, when seen from the point of vue of the pho- 



ton coupling, are precisely pole dominance and mesonic 
cloud contributions for the electromagnetic form factors. 
The objective of this paper is to set up the general for- 
malism for the evaluation of electromagnetic form factors 
in the presence of SxSB and to use it to evaluate the u, 
d anomalous magnetic moments for various models. 

In the Pauli notation, the electromagnetic current up 
to first order in the photon momentum q u is, 



euT^u = eu 



7 M +i 



/i = jU (l + a) , fi = 



2M 
e jh 
2 M c 



q v a 



(1) 



where a stands for the anomalous part of the magnetic 
moment fi and M the particle mass. The magnetic mo- 
ment of ground state hadrons is measured experimen- 
tally. For instance we have for the proton and neutron 
a p = 1.79 and a n = —1.91. In the constituent quark 
model for light hadrons we have, 



M« = \ (^d - Mu) = -0.972^ 0p 



(2) 



and the quark magnetic moments are nearly proportional 



to the charges e u = , ea 



re, which suggests that 



the gyromagnetic factor 2(1 + a) is nearly flavor indepen- 
dent. The quantity which can be measured is M / (1 + a). 
For quark flavors u and d we have, 

M u ~ (1 + a„)338 MeV , M d ~ (1 + a d )322 MeV (3) 

The constituent quark model can be applied to fit the 
hadron spectrum, with a confining interaction, an hyper- 
fine interaction, and a zero point energy, M The required 
parameters are of the order of a s = 0.974 , M u ~ Md = 
420 MeV which would suggest a sizeable a of the order of 
.15 to .3. It is also clear that we will need ad — a u ~ .05 
in order to recover the isospin symmetry. 

The Ward Identity, 

iq^S{p + q/2)T»S(jp - 9 /2) = S(p + q/2) - S(p - q/2) 
q u r»= lS-\p + q/2) - lS-\p - q/2) (4) 

is obeyed both by the bare vertex Tq and by the Bethe- 
Salpeter vertex T M M . We will show that in the limit of 
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small momentum q, this identity has the following solu- 
tion for the vertex, 

r"(p, q) = i^S-'ip) + q v T v »(p) + o(q 2 ) (5) 

where q u T v ^(p) is defined as the kernel which is not de- 
termined by the Ward identity, 



g„ [q v T^{p)\ = . 



(6) 



The Ward identity ensures that charge conservation sur- 
vives renormalization. However it does not constrain the 
kernel, which is a signature of the renormalization. In 
particular the kernel contributes to the anomalous mag- 
netic moment of fermions. 

This can clearly be seen in QED where the infrared and 
ultraviolet divergences can be removed from the photon 
propagator, 



(7) 



(pi _ p)2 {pi _ p^2 _ \2 ( p l _ p^2 _ j\2 

The vertex is given by, 

= T£ - id^T, + q v 7^ (8) 

and , up to first order in a, the contributions from the 
self energy and the kernel to the anomalous magnetic 
moment are respectively, 







)- 




) 2ir ' 





In the case of actual QED, where A — > 0, A — > oo, they 
are both infrared divergent but their sum is finite: a/2ir. 

As for QCD, there has been a considerable effort on 
how to derive quark models by integrating out, under 
various approximations, the gluonic degrees of freedom. 
An interesting and promising approach is provided by 
the cummulant expansion of the interaction term of the 
QCD Lagrangian ||. A non-local Nambu Jona-Lasinio 
type Lagrangian (NJL) is obtained when we retain only 
bilocal correlators. Therefore we hold the view that such 
quark models are appropriate to study electromagnetic 
properties of hadrons, even for light quarks, provided we 
have small enough photon momenta and the physics of 
chiral symmetry breaking is treated correctly. Therefore, 
at this stage, rather than focusing on a specific example 
of NJL we will study the static electromagnetic properties 
of a wide class of quark effective quartic interactions. 

In quark models with dynamical SxSB, the vector ver- 
tex T M is a solution of the Bethe-Salpeter equation, 



in 

P2 

Pl=P- 




where the strong interaction, which is described by a dot- 
ted line in the diagrams is iterated to all orders in the 
Bethe Salpeter equation. This equation can be written, 



v(p, q ) =rff-i 



d 4 p' 



■v(p'-p,p' + p : q) n„s(pi) 



rV, §)S(p£)n„ - v (q, P ' + P , -p' + p) n a 



tr{s( P ' 1 )r^p',q)s( P ' 2 )n a } 



(ii) 



where the —1 factor from the fermion loop was included 
in the tadpole term. The momentum dependence of the 
potential is only assumed to conserve the total momen- 
tum, and in this case it depends on 3 momenta. The 
Dirac, flavor and color structure of the interaction is de- 
termined by the f2 a matrices. In order to have dynami- 
cal SxSB, we require this structure to be chiral invari- 
ant. Substituting the Ward Identity in the ladder Bethe 
Salpeter equation for the vertex we get, 

iS^ipi) - iS^ 1 {p 2 ) = iS 1 {pi) - iS 1 (p 2 ) 
d 4 p' 



(2n> 



-v{p -p, P ' + p, q)n a [5 (pi) - s(p 2 )] n a 



-V (q,p'+p, ~p'+p) n a tr{[S(p[) - S(p' 2 )} Q a } . (12) 

For particular cases of the potential V(p[ —p\ , p' x +p 2 , p[ — 
p' 2 ) we recover the BCS mass gap equation, 



p' 



p-p' 



(13) 



provided that either the rainbow diagram vanishes or, 

V(p' ~p,p'+ p, q) = V(p'- p, p' +p,0), (14) 
and that either the tadpole diagram vanishes or, 

V(q,p'+p,-p'+p)=V(Q,p'+p,-p'+p). (15) 
Equation (|TJ) can be written, 



iS-\ P ) =iSv\p)- 



(27T)" 



-v{p' - P ,p' + p 7 o)n a s{ P ')n a 



-v (o, P ' + P , - P ' + P ) n a tr {s( P ')n a } . 



(16) 



Now we insert the expression d) for T' 1 in the Bethe 
Salpeter equation (^), in order to find the kernel qT 
and expand it up to first order in q. The equation for the 
tensor T, which is antisymmetric, is then, 



d 4 p 

(2tt) 4 



vn a (i-tr){ST^sn a } 



d 4 p 



Vn a {l-tr){J v ^ a } 



2 J (2tt) 4 

J v » = d v {S) S^d^iS) - d»{S) S^d^S) (17) 
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This is a self consistent forced linear integral equation. 
Let us consider a general quark propagator, solution of 
the mass gap equation, of the form, 



iF(p) 



t> - M{p) 

where p—^Jp^p^. The integrand J v ^ is then, 



(18) 



J"* = -i 



F 



(pa _ m 2 ) 2 



p 



(19) 



where the dot superscript denotes d/dp. In general we 
find, 

= t 1 (p){6tf,<f]} + t 2 (p) M[<f, 7 1 

+t 3 (pwiAr}-p^Y}) (20) 

Up to o(q 2 ) the electromagnetic current of the quark is 
then, 



J F 



e(l - M) 



F 



— M — (J>— M) — — + jyT"" 
p p b 

T fiU 



Y l + a [i 



2M 



M + AM 2 F(2t 1 + £ 2 ) 
1 - M 



(21) 



where the mass shell condition p = M was used together 
with the Gordon identities. The anomalous magnetic mo- 
ment a turns out to be independent of £3 and F. However 
the dependence on M is crucial in models where £1 and £2 
are finite. In those models a can be thought as a measure 
of SxSB. The quark condensate (qq) is also a functional 
of the dynamically generated mass, 



(93) 



Tlr tr 



d 4 P 

W) 4 



S(p) 



(22) 



where the trace sums colors with n c = 3, but the flavor 
is kept fixed. Thus, at the onset of the spontaneous xSB, 
we will obtain an implicit relation between a , (qq) and 
the constituent quark mass, which were simultaneously 
vanishing before the occurrence of this phase transition 
and now become non-zero. 

We will now compute F, M, a and (qq) in particular 
models which are paradigmatic cases of chiral symmetry 
breaking and comply with the constraints of the Ward 
identity. 

Model I is the first original NJL model Q. The La- 
grangian of model I is, 



qi flq + G (qq) - (qjsq)' 



(23) 



where Ci is specific to the case of 1 flavor, but its re- 
sults are similar to the ones of flavor symmetric [/^(n/) 
extended NJL models. The equations will be solved in 
the momentum representation. As usual the integrals 
are done in Euclidean space. A momentum cutoff A is 
included in order that the integral in the loop momen- 
tum is finite. Since the cutoff cannot be adscribed to the 
potential which has to be constant in momentum space, 
is must be included in the propagator, 



S( P ) 



F(p) 



JP-M + ; 



F(p) -> QEucUMan^ ~ p) ■ (24) 



With a constant potential and this momentum cutoff, 
the loops turn out to be constant, independent of the 
external momentum p. It is convenient to evaluate the 
integrals, 



d 4 p 



iF 



16n 2 



h = 



(2tt) 4 (p 2 - M 2 ) 

A 2 - M 2 ln ^1 

d A p F 
(2tt) 4 (p2 _ M 2 ) 2 



A2_ 

M 2 



16tt 2 



A 2 



A 2 + M 2 



In 1 



A*. 

M 2 



(25) 



where the solid angle 2-7T 2 is included. The mass gap 
equation is, 



p - M = i> - 2G 



d 4 p 



iF 



(2tt) 4 (p 2 - M 2 ) 
j 5 (4) + M)y 5 -tr{fi + M} 



(fi+M) 



(26) 



With the solutions M = or l = 8n c Gh(M, A). The 
parameters A and G are determined once the quark dy- 
namical mass and the quark condensate are fixed. We 
now study the kernel in model I. Because the integrals 
are constant, the antisymmetric tensor T is independent 
of p. Thus has to be of the £2 type, proportional to 
\~f 1 7 M ]- Including the structure factors f2 a we find that 
the tadpole-like term vanishes since a vil and a vlL ^ have 
a null trace. In this case of model I the rainbow diagram 
also cancels since the structure 1 ® 1 — 75 ® 75 projects 
on the terms with an odd number of Dirac 7 matrices, of 
type £1 but \~f , 7 M ] is even. Thus model I produces no 
kernel for the vector vertex and no anomalous magnetic 
moment for the quark Q. 

Model 77 is the second original NJL model ||. The 
Lagrangian is, 



C n = qi ftq + G (qq) - (q^rq) 



(27) 



where Cu is used for 2 flavors u and d. It only has an 
SU(2)a symmetry and breaks U(1)a from the onset. Its 
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results are similar to those of flavor symmetric SUA(rif) 
extended NJL models. The anzats for the propagator is 
that of Eq. ( |l8| ) , and model /J only differs from model I 
in the algebra. The mass gap equation is changed since 
t ■ f = 3 in the fermion line. We get, 



M 



2G 



d 4 p 



iF 



(2tt) 4 (p 2 - M 2 ) 



(£> + M) 



-75(i> + M) 7 5-*r{i> + A0 (28) 



=> M = or 1 = 2 2- 



l + 4n/n c Gh(M,A) 



were n/ and n c stand respectively for the number of 
flavors and colors. The rainbow diagram contributes 
in this case. The tadpole diagram contribution also 
changes. The preferred values for the parameters A and 
G are A = 1.65GeV and G = 1.23GeV~ 2 which yield 
M = MGeV, (qq) = -(.25GeV) 3 and f„ = .09GeV. 
As in model /, the tadpole will not contribute to the an- 
tisymmetric tensor which will be again of the ti type, 
T 1 ^ = t% [y u , 7 M ] . The first order term is a function of, 



d 4 p 



J 



-h M [y,y] 



(29) 



In order to evaluate the higher order terms, we calculate, 
d 4 p 



(27T)' 



-S M y] S = -iM 2 I 2 q v M f] (30) 



In this case we have two flavors with two different charges 
e u = | , e d = — \ , and two anomalous magnetic mo- 
ments af, 



e ut u — G I2 
(u «-> d) 



V 2 



M 2 e„L 



M)(f 



M 2 e d t d 



The natural parameter is 2G M 2 J 2 (A,M) = 0.004. in- 
verting this equation we find the solution, 



a u ~ -2(2GM 2 I 2 )— = 0.004 

67/ 



ad 



-2{2GM 2 I 2 ) 



0.016 



M„ = 339 MeF 
M d = 327 Af e y (32) 



Although this effect is small, it has the right sign to cor- 
rect the M u and M d inversion. If the tadpole term was 
removed from the mass gap equation then the a d — a u 
would be bigger. This is possible for instance when the 
potential has a A. A dependence, being A the Gell-Mann 
matrices. 

Model III is the simplest QCD inspired model. The 
Lagrangian is, 

1 A f A 

Cm = qi ftq+-q(x)j a -q(x)- d i yV(x-y)q(y)j a -q(y) 



In the case of model 177, the Dirac structure 7^ ® j u 
is U~A{nf) chiral invariant. For V(p) we will choose a 
color confining square well potential because of its calcu- 
lational simplicity. 

V(P - P) = ~G OEuclidian(k - \p' - p\) , (33) 

The mass gap equation is, 

t>-M ^ 4i f d*p' Sri -2j/ + 4M 
^V(p~p)F— 2 



(2tt)< 



p' 2 - M 2 



which includes the color factor of |. We now calculate the 
kernel. The first order term for the kernel is a functional 
of, 



(p 2 - M 2 y 



(34) 



The terms with two gamma matrices, of the form are 
now cancelled by the 7 a <g)7 Q of the interaction, and only 
the t\ type term remains. Therefore this model differs 
from the previous ones insofar it covers the form factor 
t\. The self consistent equation for the antisymmetric 
tensor T will also close, 

7 a s { a tf, y]} s la = 2F 2 p + {i>, [y , 7 ^]} 

\p l — M z ) 

We get, 

*) = '<--?/^r^'-")^V) 



p' 2 + M 2 



(31) to {p) = --% 



(p' 2 - M 2 ) 



rt(p') > 



-V(p'-p) 



p'-p F{p>) 



3./ (2tt)*' ^ " p 2 y 2 -M 2 ) 2 



(35) 



For the euclidean integration it is convenient to evaluate 
the angular integrals, 

I 3 (p',p) = gJ dw9 ^A -\]p' 2 +p 2 - 2wp'p^j 

= G(l + J 6 )0(1 - h)9{l + h) + 2G9{I 6 - 1) > 
h(p' ,p) = G J dwwe(k- \jp' 2 +p 2 - 2wp'pj 

[l-h)6{l+h) > 

k = '—^T, ( 36 ) 

2p'p 

and the nonlinear integral mass gap equation for F and 
M, the integral for t 0l the linear integral equation for T, 
and the integral for (qq) can be solved simultaneously, 



1 - T 2 



A 2 - p' 2 - p 2 



4 



F(p)= I 

M(p) = F(p) I dp 



dp , W,P) p' — m 



o r 6tt 2 p p' 2 + M 2 (p') 
,W,p) p' 3 F(p>) 



3tt 2 p' 2 + M 2 (p') 



-M(p') 



to(p) 



dp 



,I 5 (p',p) p' 4 F(p') 



o ' [p> 2 + M 2 {p')] 2 



t(p) = t - dp' 



m = 



dp 



6ir 2 p [p' 2 + M 2 (p'))' 

3 p 3 F(p)M(p) 



2tt 2 p 2 + M 2 (p) 



t(p') 
(37) 



The mass term has a trivial solution M = and an- 
other solution which breaks spontaneously chiral sym- 
metry. A dimensional simplification occurs if we work 
in units of A = 1. In this case the only parameter is 
G which is now adimensional. We find a critical value 
G c = 132 above which chiral symmetry occurs. In Fig. 
1 we depict the values of M, (qq) and a. We solve the 
integral equations numerically for F , M and t with the 
Gauss iterative method and using the Gauss integration 
p0| . We find that at p 2 = —1 these functions decrease 
by a factor of just .9 — > .7. Since we cannot continue 
analytically the numerical solution we use the approxi- 
mation of nearly constant F, M and t and compute the 
mass and the anomalous magnetic moment for p = 0. 
The literature prefers a < qq >= —(0.25 GeV) 3 . A dy- 
namical quark mass M — 0.33 GeV would correspond to 
G = 245A~ 2 , A = 0.74 GeV and a = 0.15 . If we now 
consider a M = (1 + a) 0.33 GeV then the lowest possible 
condensate is < qq >= — (0.28GeF) 3 which corresponds 
to G = 300A~ 2 , A = 0.69 GeV, M = 0.42 GeV and 
a = 0.28 , see Fig. 1. 

NJL models I and II are the simplest models with chi- 
ral symmetry breaking. In the NJL model I the anoma- 
lous magnetic moment a vanishes. In the model II the 
C/(l) breaking interaction yields a too small a, which nev- 
ertheless provides an example of an isospin dependence 
for a and, therefore, contributes to the u — d mass 
inversion. The reason for the smallness of the anoma- 
lous magnetic moment stems from the presence of tad- 
pole contributions and were not for this contribution and 
we would have obtained a much larger a. This is precisely 
the case of model III where a larger a is derived, com- 
patible with the nonrelativistic constituent quark models. 
We also find that M, a, < qq >, are functions of (G — G c ) 
with critical exponents which are respectively 1, 2 and 
1. The present work constitutes a first step on a more 
elaborate model unifying hadronic spectroscopy (includ- 
ing decay widths) and the electromagnetic form factors 
which have been shown to be consistent with the simple 
quark constituent picture precisely because of S\SB. 
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FIG. 1. Functions of the quark dynamical mass M, anoma- 
lous magnetic moment a and quark condensate as functions 
of the adimensional coupling G for model III 
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